The Sub-Luminal Parsec-scale Jet of Mkn 501 
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ABSTRACT 

We have combined Very Long Baseline Interferometry (VLBI) data from several programs 
in order to resolve differences in reported parsec-scale jet speeds for the TeV gamma-ray source 
Mkn 501. Data from the VLBA 2cm survey, and 8 and 15 GHz data from the Radio Reference 
Frame Image Database, have been combined with data from a 5 GHz VSOP Space VLBI observa- 
tion to determine the apparent motions of jet components in this source. The combined data set 
consists of 12 observations between April 1995 and July 1999. Four jet components are detected 
at most epochs, all of which are clearly sub- luminal (i.e. with apparent speeds < c), and two of 
which appear stationary. The established TeV gamma-ray sources Mkn 501 and Mkn 421 thus 
both have sub-luminal parsec-scale jets, in contrast to the apparently super-luminal jets of the 
majority of >100 MeV sources detected by EGRET. No new VLBI component has emerged from 
the core following the extended TeV high-state in 1997, again in contrast to the general behavior 
of >100 MeV gamma-ray sources. 

Subject headings: BL Lacertae objects: individual (Mkn 501), techniques: interferometric 



1. Introduction 

The BL Lac object Mkn 501 (1652-1-398, J1653-f 3945, 
DA426) has been well-studied at radio wave- 
lengths (e.g. Mufson et al. 1984; Gabuzda et 
al. 1992), but interest in the source was rejuve- 
nated by the discovery of TeV gamma-ray emis- 
sion (Quinn et al. 1996). A prolonged high- 
state at TeV energies in 1997 included activ- 
ity on timescales of several hours, implying the 
TeV gamma-rays originate in a relatively compact 
area (see, e.g., the review of Catanese & Weekes 
1999). The detection of correlated X-ray and TeV 
gamma-ray variability in the other well-studied 
TeV source, Mkn 421, is strong evidence in fa- 
vor of the X-ray emission being the high end of 
the synchrotron component of the spectral energy 



distribution, with the TeV emission arising from 
inverse Compton scattering of photons by the syn- 
chrotron emitting electrons, and such a model is 
also widely accepted for Mkn 501. 

On the parsec scale, VLBI observations have re- 
vealed that a jet emerges from the core at a posi- 
tion angle of ~180°, and bends by ~90° within the 
first ^2 milli-arcseconds (mas). The jet extends 
to the east until ~20mas from the core, when it 
bends further, finally reaching the position angle 
of ^^45° seen on the kilo-parsec scale (Conway & 
Wrobel 1995; Giovannini et al. 1999). The parsec- 
scale jet of Mkn 501 is one-sided, and it is assumed 
that this jet is relativistically doppler boosted, 
while the counter-jet is doppler de-boosted to such 
an extent it is invisible at the current sensitivity of 
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VLBI observations. The jets are believed to orig- 
inate in the accretion disk surrounding a central 
super-massive black hole, which for MknSOl has 
been suggested to have a mass of lO^-^siO-Si jyj^ 
(Falomo, Kotilainen & Treves 2002). 

A number of different motions have been re- 
ported for components in the Mkn501 jet, ranging 
from 0.27±0.02masyr-i (Gabuzda et al. 1994) to 
2.4masyr~^ (Giovannini et al. 1999). These stud- 
ies have relied on only a small number of epochs, 
typically less than four, and, in hindsight, have 
probably under-estimated the errors in locating 
jet components. As illustrated by Finer et al. 
(1999) for Mkn421, reliable determination of com- 
ponent motions generally requires larger, multi- 
epoch data sets. 

Here we study the parsec-scale jet of MknSOl 
from twelve VLBI observations spanning 4.28 
years, at frequencies of 5, 8 and 15 GHz. MknSOl 
lies at a redshift of 0.034 (Wills & Wills 1974) 
which, for the value of Hq = 6Skms~^Mpc~^ 
adopted throughout his paper, corresponds to a 
distance of ISSMpc. At this distance, an angu- 
lar separation of 1 mas corresponds to a projected 
linear distance of 0.72 pc. 

2. VLBI observations 

We have compiled data from three programs for 
this study: the Very Long Baseline Array (VLBA) 
2 cm survey Kellermann et al. (1998, see also 
http://www.cv.nrao.edu/2cmsurvey), the Radio 
Reference Frame Image Database (Fey & Chariot 
1997, see also http;//rorf.usno. navy.mil/rrfid.shtml), 
and a single VLBI Space Observatory Programme 
(VSOP; Hirabayashi et al. 1998, 2000) observa- 
tion. The data used in this study are summarized 
in Table 1. 

The VLBA 2 cm survey is being undertaken 
at multiple epochs to study the properties and 
evolution of over 100 active galactic nuclei. The 
MknSOl observations consisted of eight scans at 
one hour intervals of typically S minutes duration. 
Data were recorded with a bandwidth of 64 MHz 
using 1 bit samples and left-circular polarization. 
An image derived from the March 1997 data used 
in this paper was published by Kellermann et al. 
(1998). The full ten-station array was used at all 
epochs except the last, for which the Saint Croix 
telescope was unavailable. 



The Radio Reference Frame Image Database 
(RRFID) of the U.S. Naval Observatory (USNO) 
is a program to regularly image the radio sources 
used for precise astrometry. The VLBA was used 
at all epochs, although the October 199S observa- 
tion was made without the Mauna Kea and North 
Liberty telescopes. The June 1998 observation 
was made with the addition of the Fairbanks 26 m 
(Alaska), Green Bank 20m (West Virginia), Ko- 
kee Park 20m (Hawaii), Medicina 32m (Italy), 
Ny Alesund 20m (Norway), Onsala 20m (Swe- 
den), and Westford 18 m (Massachusetts) tele- 
scopes. Typically, four scans of ^3 minutes were 
made, with bandwidths of 16 MHz for the first two 
epochs, 32 MHz for the 8 GHz observations of the 
last two epochs, and 64 MHz for the 15 GHz obser- 
vation in April 1996. Right circular polarization 
is recorded for all RRFID observations. An im- 
age from the 1995 April 12 epoch at 8.4 GHz was 
published by Fey & Chariot (1997). 

The SGHz VSOP observation, in April 1998, 
was made over a 13 hour period with the HALCA 
satellite, the VLBA and the Effelsberg 100 m (Ger- 
many) telescope. Interferometric fringes to the 
satellite were detected from tracking passes to- 
taling 7 hours. In the standard VSOP observing 
mode, 32 MHz of two-bit sampled, left circular po- 
larization data is recorded. Although made at the 
lowest frequency considered here, the long base- 
lines to the orbiting telescope result in the syn- 
thesized beam-size for this observation being the 
smallest of these data. The VSOP data considered 
here were combined with a 1.6 GHz VSOP obser- 
vation one day later to derive a spectral index map 
of the source (Edwards et al. 2000a). 

3. Analysis 

The data were fringe-fit in AIPS and we have 
imaged all data ourselves using the Difmap pack- 
age. Beams were calculated using natural weight- 
ing (uvweight=0,— 1 in Difmap), with the excep- 
tion of the VSOP observation for which uniform 
weighting (uvweight=2,0 in Difmap) is more ap- 
propriate (Hirabayashi et al. 2000). An image 
from the 15 GHz RRFID observation in April 1996 
is shown in Figure 1. 

Model-fitting of the images was carried out in 
Difmap. As inspection of Figure 1 reveals, at most 
epochs four jet components were required in addi- 
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Table 1: Observations of Mkn 501 



Frequency Synthesized Beam'' Image r.m.s. Noise 



Date 


(GHz) 


(mas) 


(mas) 


(°) 


(/iJy beam ^) 


Program 


1995 Apr 8 


15 


0.96 


0.52 


-5.2 


390 


2 cm survey 


1995 Apr 12 


8 


1.33 


1.26 


18.4 


675 


RRFID 


1995 Oct 17 


15 


1.24 


0.59 


-1.6 


1195 


RRFID 


1995 Dec 15 


15 


0.96 


0.54 


-7.1 


320 


2 cm survey 


1996 Apr 23 


15 


0.85 


0.66 


2.2 


565 


RRFID 


1996 Apr 24 


8 


1.89 


1.22 


19.9 


445 


RRFID 


1996 Jul 10 


15 


0.96 


0.53 


-1.4 


275 


2 cm survey 


1997 Mar 13 


15 


0.95 


0.54 


-8.8 


205 


2 cm survey 


1998 Apr 7 


5 


0.58 


0.23 


21.9 


570 


VSOP 


1998 Jun 24 


8 


0.94 


0.76 


-19.7 


515 


RRFID 


1998 Oct 30 


15 


0.90 


0.53 


2.9 


300 


2 cm survey 


1999 Jul 19 


15 


1.26 


0.61 


-8.7 


290 


2 cm survey 



"Numbers given for the beam are the FWHMs of the major and minor axes in mas, and the position angle of the major axis in 
degrees. 



tion to the core to provide a good representation 
of the data. We have labeled these CI through 
C4, with CI being the component farthest from 
the core. Circular gaussian components were fit- 
ted at all epochs. Full details of the model-fits are 
given in Table 2. Reduced values for the fits 
are not given in the Table, as they are dependent 
on the way the data from the different programs 
was reduced. Thus, while model fits represent a 
minimum in for the given number of compo- 
nents, the comparison of x^ values between data 
from different programs is potentially misleading 
(see also Finer et al. 1999). 

The positions of all model-fit jet components 
arc plotted as a function of time in Figure 2. In or- 
der to determine the uncertainty in the component 
location for the purposes of determining compo- 
nent motions, we have projected the beam major- 
axis onto the line joining the component and core 
and then taken a fraction of this projected length 
as the error in position. For the extended, outer- 
most component, CI, we conservatively adopted 
half a projected beamwidth for the uncertainty. 
For C2 and C3 we used one quarter of the pro- 
jected beamwidth, and for the innermost compo- 
nent, C4, we used one-eighth of the beamwidth. 
Motions were determined by weighted linear fits to 
the data. As shown in Figure 2, both CI and C4 



show little evidence of motion over the 4.28 year 
period. Both C2 and C3 show clear evidence of 
motion, with an apparent component speeds for 
C2 of 0.6±0.1c, and for C3 of 0.3±0.1c. The ap- 
parent speeds along the jet, as opposed to radial 
separations, are, within errors, the same. These 
speeds supersede the preliminary values reported 
in Finer et al. (2002). 

4. Discussion 

Our component locations agree well with those 
reported from a contemporaneous 5 GHz observa- 
tion made in June 1996 as part of the VLBA Prc- 
launch Survey (Fomalont et al. 2000). We can 
also extrapolate our derived motions and compare 
them with model-fits to the 5 GHz observations at 
epochs 1987.4 (Gabuzda et al. 1992), and 1989.3 
(Gabuzda et al. 1994). The extrapolated motion 
of C3 is consistent with the positions of the K2 of 
Gabuzda et al. (1992) and Gabuzda et al. (1994), 
assuming uncertainties -^3 times larger than the 
±0.1 mas adopted by these workers. The Kl of 
Gabuzda et al. (1992) and Gabuzda et al. (1994), 
with similarly increased uncertainties, is consis- 
tent with the extrapolated motion of C2, particu- 
larly if the speed of C2 lies at the lower end of the 
range determined in §3 (assuming a constant mo- 
tion). These identifications were qualitatively sug- 
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Table 2: Gaussian Model Fits to Source Components 





Frequency 


Component 






PA" 


a" 


Epoch 


(GHz) 


ID 


(mJy) 


(mas) 


(deg) 


(mas) 


1995 Apr 8 


15 


Core 


489 






0.17 






C4 


114 


0.75 


172.1 


0.63 






C3 


84 


2.32 


147.7 


1.25 






C2 


57 


4.01 


133.1 


1.95 






CI 


(i2 


7.15 


113.0 


2.58 



"Flux density in milliJanskys. 

''r and PA are the polar coordinates of the center of the Gaussian relative to the core. Position Angle is measured from north 
through east. 

■^a is the FWHM of the Gaussian. 
[The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 



gested by Edwards et al. (2000b), but are quanti- 
tatively borne out by the fuller analysis presented 
here. 

The speeds of C4 and CI are formally consis- 
tent with zero, i.e. they appear to be stationary 
components. Before considering this further, we 
reconcile this result for C4 with the observations 
of Marscher (1999), who reported the detection 
at 22 GHz of a resolved component between 0.5 
and 1 mas from the core, with an apparent mo- 
tion derived from three epochs between April and 
August 1997 of 0.96±0.1 masyr"'^, corresponding 
to 2.3±0.2c. This location is consistent with our 
C4, however we do not see such rapid motion over 
the four year period. The motion reported by 
Marscher corresponds to 0.27 mas in the 0.29 years 
the observations spanned. In our data, C4 ranges 
between 0.59 mas from the core (1998 April 7) and 
0.83 mas from the core (1999 July 19), a range of 
0.21 mas, similar in magnitude to that of Marscher 
(1999). Any attempt at further interpretation is 
complicated by the fact that there are likely to 
be frequency-dependent offsets in the separation 
of components from the core (see, e.g., Lobanov 
1998) in our data, which would be most impor- 
tant for C4. 

Stationary components have been reported for 
a number of sources in the past, with a detailed 
study being made as part of the multi-epoch mon- 
itoring program of Jorstad ct al. (2001a). The 
monitoring revealed that the super-luminal speeds 
detected for these of EGRET-detected blazars 
were much faster than for the general population 



of bright compact radio sources, however evidence 
was also found for at least one stationary com- 
ponent in 27 of the 42 sources (Jorstad et al. 
2001a). Jorstad et al. suggested that the station- 
ary components within several parsecs of the core 
were associated with standing recoUimation shocks 
caused by pressure imbalances at the boundary 
between the jet and the surrounding medium. In 
contrast, the stationary components further from 
the core tended to be associated with bends in 
the parsec-scale jet. There is support for this sce- 
nario in our data. C4 is located at a projected dis- 
tance of ~0.5 pc from the core, and is quite plausi- 
bly associated with a recoUimation shock. CI, on 
the other hand, is an extended component, which 
1.6 GHz VLBI imaging has revealed is associated 
with a significant change in the jet from a bright 
"spine" to a limb-brightened morphology (Giovan- 
nini et al. 1999). 

If we assume our fastest observed pattern speed 
(0.6c for C2) reflects the bulk apparent speed of 
the jet, then we can solve for the intrinsic speed 
and angle to the line-of-sight, provided we also 
have an estimate of the Dopplcr beaming factor. 
A Dopplcr factor 5 ^ 10 is inferred from the 
TeV observations of this source (e.g., Tavecchio, 
Maraschi, & GhiseUini 1998). The VSOP observa- 
tions yield our best measurement of the radio-core 
brightness temperature, 4xlO^^K. This is consis- 
tent with a Dopplcr factor of ^10 if the source 
is in equipartition (Readhead 1994), but it is also 
consistent with lower Doppler factors if equiparti- 
tion is violated (Kellermann 2002). If we accept 
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the values of 0.6c and 10 for the apparent bulk 
speed and Doppler factor, then the Lorentz factor 
of the MknSOl jet is 7 = 5 (?; = 0.98c) and its 
angle to the line-of-sight is 6 = 0.7°. Such a small 
angle to the line-of-sight may be expected of a 7- 
ray blazar, although subluminal apparent speeds 
are in general not expected (see the Monte Carlo 
simulations of Lister 1999). 

Alternative kinematics that do not place such 
tight constraints on the angle to the line-of-sight 
assign the Doppler factor measurement to the 
TeV-emitting region (on the light-day size scale) 
and the apparent bulk speed to the VLBI jet (on 
the light-year size scale), and allow a change in 
the bulk Lorentz factor or angle to the line-of-sight 
in the intermediate region. If the jet in the TeV- 
emitting region has, e.g., ^ = 5° and 7 = 7 (enforc- 
ing 6 = 10), then a decrease in the Lorentz factor 
to 7 = 2 would reproduce the observed appar- 
ent speed in the VLBI jet. Such a deceleration of 
electron-positron jets close to the core is proposed 
by Marscher (1999) for the TeV blazars. A change 
in angle to the line-of-sight, perhaps accompany- 
ing the large bend in the jet seen ~2 mas from the 
core, cannot by itself reproduce the observed val- 
ues; a jet with S = 10 has a minimum Lorentz 
factor of 5, and a jet with 7 = 5 can only have an 
apparent speed of 0.6c in the large-angle solution 
for 9 > 90°. Any set of kinematic parameters must 
also be constrained by the one-sided appearance 
of the source; the example above with 9 = 5° and 
7 = 2 would have a jet-to-counterjet brightness ra- 
tio greater than ~200, somewhat higher than the 
limit that can be placed from our observations. 

Similar values for the Doppler factor and ap- 
parent jet speed apply to the other well-studied 
TeV blazar, Mkn 421 (Finer et al. 1999). Prom 
these two sources, it appears that TeV blazars as 
a class may either have very small angles to the 
line-of-sight {9 < 1°), or may decelerate signifi- 
cantly between the TeV-emitting region and the 
parsec scale. 

It is notable that no new component has 
emerged from the core after the prolonged TeV 
high-state in 1997. A component with a speed 
similar to that of C2 or C3 would now be ~0.5 mas 
from the core and would have been detected at the 
latter epochs. This suggests that events which give 
rise to extended TeV (and associated X-ray) activ- 
ity are different in nature to those which result in 



the production of new VLBI components (see also 
Marscher 1999). Mkn 421 and Mkn 501, for which 
the inverse-Compton component of the spectral 
energy distribution (SED) peaks at TeV energies, 
have sub-luminal component speeds and appar- 
ently no new component emerging after epochs 
of TeV activity. In contrast, sources with the 
inverse-Compton component of the SED peaking 
at GeV energies tend to have the emergence of 
new, super-luminal, VLBI components associated 
with GeV flaring states (Jorstad et al. 2001b). The 
detection of more TeV gamma-ray sources by the 
next generation of air Cerenkov telescopes will en- 
able these apparent trends to be investigated more 
quantitatively. 
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15 GHz 



Core 
C4 



C2 



C3 



C1 



Relative Right Ascension (mas) 

Fig. 1. — Image of Mkn501 at 15 GHz from the 
RRFID observation in April 1996. The posi- 
tions of the core and jet components arc indicated 
(see Table 2 for details). The beam, 0.85 mas 
X 0.66 mas (FWHM) at a position angle of 2°, 
is shown at bottom left. The contours arc —1 
(dashed), 1, 2, 4, 8, 16, 32 and 64% of the map 
peak of 485 mJy/beam. 




M t 



■ - — i s 

1997 1998 1999 

Time (years) 



Fig. 2. — Component positions and weighted lin- 
ear fits to component motions. The squares denote 
15 GHz observations, triangles denote 8 GHz ob- 
servations, and diamonds are used for the 5 GHz 
VSOP observation. Note that the two observa- 
tions in April 1996 have been offset from each 
other in the plot for clarity. 
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